Scaling of excitons in carbon nanotubes 
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Light emission from carbon nanotubes is expected to be dominated by excitonic recombination. 
Here we calculate the properties of excitons in nanotubes embedded in a dielectric, for a wide range 
of tube radii and dielectric environments. We find that simple scaling relationships give a good 
description of the binding energy, exciton size, and oscillator strength. 
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The optical properties of carbon nanotubes have re- 
ceived increasing experimental and theoretical atten- 
tion. Optical absorption and emission spectra of car- 
bon nanotubes have been studied by a number of groups 
EL IE IE B> IE El; an d electro-optical devices have already 
appeared [21 la] ■ Initial attempts to explain the experi- 
mental observations naturally took independent-electron 
theory as their starting point. However, theoretically it 
is now clear that emission is dominated by excitonic re- 
combination H EE El EE 

A number of theoretical approaches have been used 
to describe these excitons. One approach involves vari- 
ational calculations 0,0]. While valuable, these have 
been limited to an effective-mass approximation, and do 
not address issues of spectral weight. The most accu- 
rate description is provided by an ab initio solution of 
the Bethe-Salpeter equation using GW-corrected quasi- 
particle energies EE However, it is not currently feasi- 
ble to apply this computationally intensive approach to 
a wide range of nanotube sizes or environments. 

Here we use an intermediate level of theory to provide 
a broad overview of the exciton properties. We calculate 
the excitonic properties of nanotubes embedded in di- 
electric media, for the range of tube radii and dielectric 
constants most relevant to potential applications. Wc 
find that the exciton size, binding energy, and oscillator 
strength all exhibit robust (though approximate) scaling 
relationships. The relationships obtained for the exci- 
tonic properties can be used to better understand and 
optimize the operation of nanotube opto-electronic de- 
vices. 

The proper procedure for the calculation of excitons 
has been described in detail in Ref. EE • It involves solv- 
ing the Bethe-Salpeter equation, 
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where the kernel IC k ,k' describes the interaction between 
all possible electron-hole pairs of total momentum q exc , 
and Afc is the quasiparticle energy for a non-interacting 
electron and hole with wavevector k and q exc ~ k. The 
exciton momentum q exc is equal to that of the exciting 
photon, and is hereafter approximated by q exc = 0. We 
approximate the quasiparticle energies by eigenvalues of 



the tight-binding Hamiltonian 

0G3 (* = 3.0 eV), with 
any additional self-energy corrections restricted to the 
so-called "scissors operator", in which the self-energy is 
approximated by a rigid shift of the conduction band rel- 
ative to the valence band. Since the quasi-particle band- 
structures are not well known for nanotubes of varying 
diameters and embedding media, we report only proper- 
ties that are not affected by the magnitude of this shift. 

For the optically active singlet excitons, the interaction 
has two contributions, direct {JC d ) and exchange (JC X ): 
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where the direct (exchange) term is evaluated with the 
screened (bare) Coulomb interactions EE- The un- 
screened Coulomb interaction between carbon p z orbitals 
is modelled by the Ohno potential, which realistically de- 
scribes organic polymer systems: 
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where is the distance between sites i and j, and U = 
11.3 eV is the energy cost to place two electrons on a 



single site 



0). Our results are not sensitive to the 



value of U when the size of the exciton is large. 

An ideal calculation would include the nonlocal dielec- 
tric response of both the nanotube itself and the medium 
in which it is embedded. Here, for computational simplic- 
ity, we replace this complicated response function with a 
single dielectric constant e EE- This is most accurate 
for narrow tubes, and for embedding media with large 
dielectric constants. In this regime, the exciton length 
along the tube is large relative to the tube radius, and 
most of the dielectric screening occurs in the surrounding 
medium. The screening is then well described by the di- 
electric constant e of the nanotube environment. For iso- 
lated nanotubes, or tubes in low-e media, this treatment 
is not very accurate. Fortunately, it is most accurate 
in precisely the regime of greatest technological interest. 
Modulated electro-optical devices 0, IE are most practi- 
cal for relatively narrow tubes embedded in Si02 (e ~ 4) 
or higher-e materials. 

To calculate the optical properties, we evaluate the 



2 




FIG. 1: (a) Binding energy of first optically active exciton, vs. 
e, in four semiconducting zig-zag tubes: (13,0), (19,0), (25,0), 
and (31,0). (b) Scaling of binding energy of first and second 
exciton (red dots) in semiconducting tubes with all possible 
chirality (156 tubes with d — 1.0 — 2.5 nm, e = 2—15). Here R 
and m are in a.u. The black solid line is the best fit to Eq. J7J 
for e — 4 — 15 (3% RMS error over this range), corresponding 
to a = 1.40 and A b = 24.1 eV. 



imaginary part of the dielectric function for light polar- 
ized along the nanotube axis |l7j : 
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where P is the dipole matrix element [18J. The optical 
response Eq. Q is the same as derived in the presence of 
the GW non-local potential ^t|- e 2(w) obeys a sum rule, 
where J £2<lu) cx J2k Pcvifyl ^ s a constant independent 
of the strength of the screened interaction e 2 je. 

We solve the BSE equation Q by direct diagonaliza- 
tion, choosing a k sample sufficient to converge the low- 
energy optical spectra (and a fortiori the binding ener- 
gies) . We calculate the binding energy, size, and spectral 
function for singlet excitons. The binding energy of the 
first optically active exciton vs. e is shown on Fig. 
for four zig-zag tubes with diameters d — 1.0 — 2.5 nm. 
(There is another singlet state 3-5 meV lower in energy, 
but it is optically silent by symmetry.) 

The dependence of binding energy Ebi on e in Fig. 
can be fitted well with a power law, with the exponent 
being almost independent of the tube diameter. This 
suggests a more general power law scaling, which we can 
motivate in an effective mass approximation as follows. 
Given a variation wavefunction described by a single pa- 
rameter L that scales the size along the tube axis, the 
exciton binding energy is: 
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effective mass. The second term is the potential energy, 
which depends on the exciton size via the dimensionless 
function f(L/R). Then the exciton binding energy is 
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where we approximate the function h by a power law over 
the range of interest, with empirical parameters a and 
Ah. The effective mass m depends on the tube indices 
1 1 11 ] (i.e. on radius and chirality). 

In 3D semiconductors, the potential energy is oc l/L, 
and the energy is minimized when the exciton size Ls oc 
e/m, so the binding energy scales as Eb oc m/e 2 . This 
corresponds to true scaling, with a power law a = 2 in 
Eq. (7J. In the case of nanotubes, the power-law scaling 
is only an approximation. Nevertheless, for the most im- 
portant range of tube sizes and dielectric constants, the 
behavior is rather well described by a power-law scaling 
in R and e with a single value of a. Indeed, all the binding 
energies for the first and second excitons in semiconduct- 
ing tubes, with all possible chiralities, (d = 1.0 — 1.5 nm) 
collapse onto a single curve shown in Fig. U>. Similar 
energy scaling was reported by Pedersen jlOj in a varia- 
tional effective-mass model. In the range e > 4, where 
our approach is most reliable, we obtain the best fit with 
a = 1.40. 

The second exciton that is optically active derives pri- 
marily from the second band of the nanotube. It falls 
within the continuum of the first band, and so becomes a 
resonance with a finite lifetime |l2l | . By artificially turn- 
ing off the interband coupling, we determine that this 
coupling has very little effect on the exciton energy. 

Considerable attention has been focused on the ra- 
tio EbijEbi between the binding energies of the first 
and second excitons [ill Il2j |. [Ebi is defined relative 
to the second-band quasiparticle gap. Note that the 
exciton formation energies involve also the quasiparti- 
cle bandgaps.) The scaling relation of Eq. J7J predicts 
Ebi/Ebi = (m.2 / 'm\) a ~ l , where ra-i and mi are the effec- 
tive masses of the first and second bands. In the case of 
zig-zag tubes: 
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Here the first term is the kinetic energy, and m is the 



Here Ai and A 2 are the tight-binding bandgaps; a is 
the graphene lattice constant; and for tube indeces (n,0), 
(7=1 if mod(n,3)=l and o=-l if mod(n,3)=2. 

It is common to treat the gap values in the infinite- 
radius limit, A!f = 2Af = 2ta/V3R. This is rather ac- 
curate (within 5%) for tubes with d = 1.0 — 2.5. On the 
other hand, for the same range of diameters, the effective 
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mass ratio mi/m\ varies from 3.4 to 1.3, approaching the 
infinite-radius limit mi/m\ — > 2 much more slowly than 
the gap ratio. Thus caution must be used in discussing 
available experimental data in terms of the R — > 00 limit 
In particular, for the (8,0) tube ma/mi = 0.96, 
and according to Eq. J7J the binding energies of the first 
two excitons should be very similar. Indeed, the accu- 
rate first-principles calculations by Spataru et. al. 0] 
find the binding energies of the first and second exci- 
tons (A'i and C'i in pj) to be 0.99 and 1.00 eV. Using 
e = 1.93 to best reproduce this, our calculations give 
E bl = 0.99 eV and E b2 = 1.05 eV. In contrast, the (10,0) 
tube has m^/mi = 4.14, and for the same e = 1.93 we 
find E b i/E b \ — 1.41. [The simple m a ~ 1 scaling is not 
accurate for such small e.] It is important to note that 
effective mass dependence similar to Eq. (|SJ) holds also 
for chiralities other from zig-zag tubes. Thus we expect 
exciton properties in tubes of index (m,n) to depend pri- 
marily on whether mod(n-m,3)=l or 2, independent of 
the chiral angle. 

To quantify the exciton size, we use the root-mean- 
square (RMS) distance between electron and hole, L$- 
The size L\ of the first exciton is shown in Fig. as 
a function of e, for four different tube diameters. The 
size is approximately linear in e. From Eq. JSJ), L$/R is 
expected to be a function of mR/e. Combining this with 
the observed linear dependence on e, we anticipate that 
the exciton size will obey the scaling relationship: 
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This is confirmed in Fig. [2b, which shows a linear de- 
pendence of L\/R on e/miR in all semiconducting tubes 
with d = 1.0 — 2.5 nm, for all chiralities. 

The exciton size directly affects observable quantities 
such as the exciton oscillator strength and the radiative 
lifetime. The exciton oscillator strength is proportional 
to the probab ility to find an electron and a hole at the 
same position [l9( . In 3D semiconductors this is inversely 
proportional to the exciton volume 1/L 3 . In the case 
of nanotubes the electron and hole wavefunctions are 
confined in two dimensions and therefore the oscillator 
strength should be inversely proportional to the exciton 
size L\. 

Typical optical absorption spectra are shown on 
Fig. Ot-c, calculated for a (19, 0) tube in different di- 
electric media. As s increases, the spectral function con- 
verges to the non-interacting limit. For e = 10, the spec- 
tral weight transfer to the first and second excitons, as a 
fraction of the total spectral weights for first and second 
bands in the non-interacting limit, are 71% and 55 % re- 
spectively. The second exciton resonance is more bound 
than the first, by 63 meV vs. 43 meV. The higher spectral 
weight transfer to the first exciton is due not to stronger 
binding, but rather to the smaller band gap Ai. 
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FIG. 2: (a) First exciton RMS e-h separation L\ in four zig- 
zag tubes (13,0), (19,0), (25,0), and (31,0). The solid lines 
are the best linear fits. The slope k scales approximately as 
m _1 : the product kmi (in a.u.) equals 0.167, 0.160, 0.157, 
and 0.155 for tubes with diameter 1.0, 1.5, 2.0, and 2.5 nm, 
respectively, (b) The linear scaling of Li/R with e/mR in 
semiconducting tubes of all possible chirality with d = 1.0 — 
2.5 nm, and e = 2 — 15. The solid line is the best fit to 
Eq. © to the results in the range e > 4, giving Al = 2.13 
and Bl = 0.174 (for m and R in a.u.). The RMS discrepancy 
between the fit and the full calculations is 2% for the subset 
of data having e > 4 or e/mR 2 > 9.5. 



The probability argument [l|j along with Eq. sug- 
gest the following scaling relation for the exciton oscilla- 
tor strength: 
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where Iq is the spectral weight of the first band for non- 
interacting particles. The second term is a higher order 
correction due to the band mixing and the non-constant 
band-to-band optical density. Figure |3i confirms this 
scaling. While L\ is not directly observable, the figure 
is virtually unchanged if we use the actual calculated L\ 
values instead of the scaling expression Eq. J^J for L\. 
Thus 1\ / Jo obeys rather well an explicit scaling relation- 
ship with e/Rm. 

The radiative lifetime of the excitons in carbon nan- 
otubes is a key factor for possible applications in pho- 
tonics and optoelectronics. The radiative lifetime of an 
exciton is inversely proportional to its oscillator strength 
[2fj|. In the regime of large binding (e < 3), 1\ ~ Jo- 
Then the oscillator strength per atom is almost inde- 
pendent of tube diameter and chirality, and is equal to 
f w (O.OMeV- 1 ) £ oxc . 

We emphasize however that the radiative lifetime and 
luminescence efficiency of nanotubes involve other factors 
as well. In principle, a single exciton coherently captures 
spectral weight from a macroscopic region |2l| : so the 
lifetime of an exciton actually depends on the coherence 
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FIG. 3: Absorption spectra £2 in (19,0) tube in dielectric 
environment (a) e — 00 (equivalent to no e-h interaction), 
(b) e = 10, (c) e — 4. (E s is the unknown self-energy shift.) 
Note expanded scales for dotted lines in continuum region in 
(b) and (c). The fractional spectral weight transfer to the 
first exciton is I1/I0 =0, 0.71, and 0.95 respectively. Spec- 
tra are broadened with a Gaussian width of 0.0125 eV. (d) 
The scaling of spectral weight transfer to the first exciton, 
Ii, according to Eq. 1101 . in all semiconducting tubes with 
d = 1.0 — 2.5 nm and e — 2 — 15 and all possible chiralities. 
The best fit to Eq. iJTDJ (RMS difference 3.5%) is obtained 
with Ai = 1.22 eV 2 nm 2 and Bi = 1.61. 



length in the nanotube, which in turn depends on en- 
vironment and temperature. (If the coherence length is 
sufficiently large, other lengthscales such as tube length 
or photon wavelength can become important.) Another 
important factor is that electrons and holes are rela- 
tively unlikely to form optically active excitons, because 
there are far more excitons that are optically inactive. 
These include triplet and other dipole-forbidden excitons 
at lower energy than the optically active exciton. Most 
importantly, only a tiny fraction of excitons have a total 
momentum compatible with photon emission, so phonon 
scattering plays an important role [2^ . 

In conclusion, we have calculated optical spectra of 
carbon nanotubes including the electron-hole Coulomb 
interaction by solving the Bethe-Salpeter equation Q in 
a tight- binding wavefunction basis set. We find scaling 
relations with respect to the tube radius and dielectric 
constant e, for the binding energy Eq. Q, exciton size 
Eq. ©, and oscillator strength Eq. l(TUj). Thus the ab- 
sorption and emission properties depend on the dielectric 
media in which the nanotube is placed. We find a strong 
dependence on tube index (chirality) , but only via the ef- 
fective mass. This depends strongly on whether mod(n- 
m,3)=l or 2, but is otherwise insensitive to the chiral 
angle. 
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